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K^ ' Abstract. 

^~ ' We obtain and investigate the regular eigenfunctions of simple differential operators 

tH- , x^ (F^^ /dx'^'^^, r = 1,2,... with the eigenvalues equal to one. With the help of these 

04 ' eigenfunctions we construct a non-unitary analogue of boson displacement operator 

[*-.. i which will be acting on the vacuum. In this way we generate collective quantum states 

^D [ of the Fock space which are normalized and equipped with the resolution of unity 

^"^ ' with the positive weight functions that we obtain explicitly. These states are thus 

• • . coherent states in the sense of Klauder. They span the truncated Fock space without 

.5^ I first r lowest-lying basis states: |0), |1), . . . , |r — 1). These states are squeezed, are 

j^ ' sub-Poissonian in nature and are reminiscent of photon-added states at Agarwal et al. 
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1. Introduction 

We are concerned in tliis work witli a procedure of generating collective quantum states 
spanning a prescribed infinite subset of the Fock space {|n)}5^o- The basis states \n) are 
orthonormal, {n\n') = 6n,n' and are eigenstates of the boson number operator h = a^a, 
n\n) = n\n), where a^ and a are the boson creation and annihilation operators satisfying 
[a, a^ = 1 and a|0) = 0. 

Our objective is to generate in a systematic way families of states sharing their 
main properties with the so-called standard coherent states (CS). These states are the 
CS in the sense of Klauder [1]. It means that they are: a) normalized, b) continuous 
with their label (generally a complex number z) and c) they admit the resolution of 
unity with a positive weight function. 

The central tools in the domain of standard CS is the unitary displacement operator 
©(z) = exp(za''" — z*a) whose action on the Fock vacuum |0) generates the normalized 
standard CS denoted by |2;) [2] 

\z)=V{z)\0) (1) 

= exp(-|2|V2)e"'^'|0). (2) 

Note that the state \z) is the eigenstate of the annihilation operator a: 

a\z) = z\z). (3) 

We expose now a series of arguments leading to a novel approach to construct families 
of CS. This approach lends itself to many generalizations. 

First of all let us consider quantum states similar to ([1]) and ([2]) but created with 
higher powers of boson creation operator. They may be generated by powers of a^ 
appearing in the exponential, via 

|z,p)~e^('^')nO). (4) 

However, the convergence arguments forbid the creation of states \z,p) for p > 2, as they 
cannot be normalized, contradicting the requirement a) above. Thus it is impossible 
to get the higher order CS in this manner [3]. A possible way out from this dilemma 
is to attempt to use functions other that the exponential assuring that the states are 
normalizable. 

From now on, let us employ the formal operational equivalence [1] 

d 



dx' 



1 ^ [a,a'^] = 1, (5) 



which permits one to relate d/dx and x respectively to the annihilation a and creation 
a^ operators. In this sense the displacement operator acting on vacuum, see Eq.([2]) 
is equivalent to e^^', which in turn is the eigenfunction of the operator d/dx with the 
eigenvalue z, reminiscent of Eq.Q. Following this initial idea we will seek eigenfunctions 
of operators more general than d/dx and shall attempt to construct CS associated with 
them. 
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In this paper we introduce a new generalization of the boson displacement operator 
enabling to construct in the systematic way the coherent states of photon-added type 

The plan of the paper is the following: in Sec. |2]we consider the simple differential 
operators x^ d'^~^^ / dx^~^^ , r = 1,2,... and find their eigenfunctions. Then, in Sec. |3]with 
the help of these eigenfunctions and by using the relations ([5]), we construct a new 
non-unitary analogue of displacement operator. The action of so obtained displacement 
operator on |0) yields quantum states satisfying the requirements a), b) and c) above. 
We demonstrate that these CS bear analogy to the so-called photon-added states [5]. 
In Sec. m we show that these CS satisfy certain eigenequation analogous to Eq.Q 
above. Further we define physical and statistical properties of these states, like the 
Mandel parameter, the metric factor and the squeezing. The resolution of unity and 
the associated Stieltjes moment problem are considered in Sec. [5l In Sec. Owe present 
many examples of properties of CS which were introduced in Sec. HJ calculated for 
various values of parameter r. The last Sec. [7] is devoted to discussion and conclusions. 
In Appendix we discuss the nature of the Stieltjes moment problem of Sec. [51 We 
demonstrate the non-unique character of solutions of the resolution of unity obtained in 
Sec. El 

2. Eigenfunctions of differential operators 

The starting point of our subsequent construction of quantum states is the determination 
of the eigenstates of differential operators in the form x^rf'""'"^, r = 1, 2, . . ., where d = -^, 
i.e. finding the eigenfunctions (with the eigenvalues equal to one) denoted by Eg{r,x) 
satisfying 

x'd'-+^Eg{r,x) = Eg{r,x), r = l,2,..., g=l,...,r + l. (6) 

The Eq.(|6|) is an ordinary differential equation of order r + 1 and it possesses r + 1 
independent solutions labelled by q, but only one of them is Taylor expandable at x = 0, 
see [6], and precisely this one, denoted by E{r,x), is of interest for our applications. 

The implementation of Eq.([6]) via Frobenius recursive method [7] gives the explicit 
form: 

E{r,x)=x'- oFr{[ ],[2,3,...,r + l],x), (7) 

which satisfies the following r + 1 "initial" conditions at x = 0: 

^, , dPE(r,x) , d''E(r,x) , , , ^ 

^(-,0) = 0, -^-^\^^, = 0,p=l,...,r-l, -^^\.=o = rl. (8) 

In Eq.(I7]) pFq{. . .) is a generalized hypergeometric function ^M,- Note that the Taylor 
expansion of E{r,x) at x = starts with x^. Only the lowest order case r = 1 can be 
written down in terms of standard functions: 

E(l,x) = v^/i(2v^), (9) 
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where Ii{y) is the modified Bessel function of the first kind. For r > 1 the functions 
E{r, x) are related to the so-called hyper-Bessel functions [101 IHl US]- (Note in passing 
that the irregular solutions of Eq.(l6]), Eq{r,x) with \im.x^Q Eq{r,x) = oo, are explicitly 
known only in the case r = 1: Ei{l,x) = \fxK\(f2^fx\ where K\{y) is the modified 
Bessel function of the second kind. These solutions do not enter our considerations). 

3. Generation of coherent states \z)r 

We introduce the boson creation and annihilation operators a^ and a, obeying the 
commutation relation [a, a^\ = 1 and we associate with them the Fock space of 
normalized and orthogonal vectors \n), {n\n') = Sn,n', n, n' = 0,1,..., satisfying 
the usual relations a^|n) = y/n + l\n + 1) and a^a|n) = n\n) = n\n), which is the 
eigenequation of the Hamiltonian of the harmonic oscillator. 

In order to incorporate the standard CS into a framework and notation that we 
develop later in this paragraph we note that the exponential exp(2;x) can be written as 
oFo([ ], [ ], zx) and consequently we rewrite T>{z) acting on |0) via Eq.([T]) as: 

oi^o([],[],klV2)' ^ ^ ^ 

We extend this idea to higher order operators appearing in Eqs.Q and ([7]) by defining 
the analogue of the displacement operator relevant to x^d^~^^ as proportional to E{r, za^) 
given by Eq.Q. 

To do so, for a complex number z we introduce now a continuously parametrized 
family of normalized quantum states. They are defined by using a non-unitary analogue 
Vriz) of the displacement operator which is proportional to E{r,za^), acting on the 
vacuum |0): 

\z)r ^ Vrizmr = K'^M^r)^P0^ |0) (11) 

= ^' ^Jj^''^ {a^r oFr{[ ], [2, 3, . . . , r + l],za^)\0), (12) 

where b{r) is a numerical factor equal to ni=o ^' ^^^ ■^ri\z\'^) is the normalization 
function obtained from the condition j.{z\z)r = 1: 

K{x)= lr\ Ylikiy] oF2r{[ ],[l, 2,2,. . .,r,r ,r + l],x), a; > (13) 

\ fc=0 / 2r— 2 entries 

which is a perfectly converging function for all values of x [9] . Formulae (fTTj) and (TT2I1 
are key equations of our method. 

We identify the numerical factor in front of o-^2r- in Eq. (ll3l) as p.7^(0) for reasons 
which will become clear later, see Eq.( l22|) below. Note that in Eq.l TTTl) Vr{0) = 1. The 
overlapping factor r{z\z')r is equal to Mriz* ■ z'), which is everywhere non-vanishing, 
except for a set (of measure zero) of zeros of Afr{y) in the complex plane y. Some 
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informations about the localizations of zeros of the function J^riv) can be obtained from 
the Hurwitz criterion [13] but this is beyond the scope of this work. 

The expansion of \z)r in terms of basic states \n) reads from Eg. (1121) : 

oo „ 



/ |n + r), (14) 

~^ n\{n + 1)! . . . (n + r — \)\^ {n + r)! 

from which we deduce that |z),. is spanned on the Fock space without first r lowest-lying 
states: |0), |1), . . . , |r — 1). The reason for that is explained above, see the remark after 
Eq.([8]). As such, the |-2)r-'s are reminiscent of but not directly related to the so-called 
photon-added coherent states, introduced in [5] and later studied in [Ml [T5l flGl fT71 fT8] . 
For any hermitian Hamiltonian H, in particular the harmonic oscillator Hamiltonian 
H = a'^a + 1/2, the time evolution of state \z)r = \z,t = 0)r is given by \z,t)r = 
e~'^^^\z,0)r, which immediately implies that AArd^P) is constant in time. 



4. Eigenproperties of coherent states l^;), and some expectation values 

In full analogy with many known varieties of coherent states [T71 [T9| [20l [2T] the state 
\z)r turns out to be an eigenstate of a certain generalized boson annihilation operator. 
In fact we show that this operator is equal to {a^Y a''"'"^: 

{a^Ya^ ■ a\z)r = z\z)^ 



Y^a{r,k){n-\Y 



k=\ 
r+1 



\z)r = aTT(n — h)\z)i 



fc=i 



y \a{r -|- 1, A;)|(n — r — 1) 



fe-i 



fc=i 



\z\ 



(15) 
(16) 

(17) 



where cr(r, k) are the conventional Stirling numbers of the first kind [22], and n = a^a. 
In order to prove the eigenproperty of Eq. (IT5|) it suffices to observe that in the expansion 
of Eq. (1141) the use of relation 



(a^Ya'^^ln + r) 



(n 



1 + r)! 



y/n + r |n + r — 1), 



(n-1)! 

and a standard change of summation index immediately gives the required result. 
Furthermore, the Eq. flT6|) is a consequence of the generating function for a{r,k) [22] 
whereas its alternative form in Eq.l lTTl) comes from a slight modification of Eq.(ll) in 
Ref. % 

Since Eqs. (TT6|) and (TT71) imply that 

afr{n)\z)r = z\z)r, (19) 

with fr{x) = X]fc=i '^('"' A;)(x — l)'^, which is a simple function of x, the states \z)r fall 
into the category of non-linear coherent states, frequently discussed in the literature 
[201 Eal [23 [251 [26] in different contexts. 
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The knowledge of Fock-space expansion of \z)r of Eq. fll4p allows one to obtain 
explicitly the expectation values of many operators involving a's and a'^^s. Below we 
quote some of them: 



and 



(2|(aT)PaP|z) 



M{ci^Yci'\z)r 



p-r Jp 



Mrix) dxP 



X = \z\ 



p = 0,l,.... (20) 



Afr{\z\') 



n=0 



[n + r)\{n + r + p — s)\\ 
Pr{n)pr{n + p - s) 



\2n 



[U + r — S]\ 



(21) 



where p, s = 0, 1, . . . and 



pr{n) 



'r-l 



n(^+^)! 



.fc=0 



'r-1 



[n 



+ r)\, Pr(0)= n^' ^'' '^ = 1'2, 



(22) 



vfc=0 



Armed with Eqs. (l20l) and fl2T|) the following physical and statistical characteristics of 

the state \z)r can be readily calculated: 

- The probability Pr{k,x) of finding the vector \k) in the state \z)r, k > r: 



±y [K^Xj 



X 



k—r 



J\fr{x) Pr{k — r) 

The average number nr{x) of bosons in the state \z)r'- 



nr{x) = r{z\n\z)r = / ^ k Pr{k, 



x) 



(23) 



(24) 



k=r 



More generally nPr{x) = r{z\n'''\z)r, p = 1,2, . . . . 

- The Mandel parameter QM,r{x), which yields information about the deviation of the 

probability distribution Pr{k,x) from the Poisson distribution, is defined as (see Ref. 

I): 



n\{x) - K(x)]^ 



X 



nrix) 



{x'K{x))" {x'-K{x))' 



^Mr{x))' 



[X 



X^'Nrix) 



(25) 



As it is the well-known [27] in the Poissonian case, relevant for standard CS, we have 
QM,r{.x) = while for QM,r{x) < (resp. QM,r{x) > 0) we say that the distribution is 
sub- Poissonian (resp. super- Poissonian) . 
- The metric factor Urix): 



Ur{X) 



m 



X) 



X 



Mr 



X) 



(26) 



It describes the geometric features of coherent states: uJr{x) = 1 for standard CS and a 
deviation from it measures the non-standard behaviour of states in question |28j . 
- Squeezing in both X and P, where X- and P-quadratures are given in terms of a and 
a^ by usual formulae 



X = (a + a'^ 



P 



-lia — a' 



(27) 
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The uncertainties in X and P are given in terms of standard formulae [H] 

= ^ (l + 2{a^a) + (a^) + (at') - (a)^ - {a^f - 2{a){a^)) , (28) 

(AP)2=752_p2 

= ^ (l + 2{a^a) - (a^) - (a^') + (a)^ + {a^)^ - 2(a)(at)) , (29) 

where in Eqs.(l28l) and ( J29l) all the averages (■ ■ ■) are understood to be calculated in 
l^;)^. The state \z)r is called the squeezed state if the uncertainty at least in one of the 
observables X or P is less than 1/2 |29j . 



5. Resolution of unity and the Stieltjes moment problem 

Since the first r Fock states are absent from l^)^, the appropriate unity operator relevant 
for this situation is [151 [181: 



(30) 



Ir = y^ |^)(^| = /, liT- + r){n + r\ 



n=0 



which is an infinite sum of orthogonal projection operators, and the resolution of unity 
for the states \z)r will read 



» oo 

/ d'^z \z)r Wrilzl"^) r{z\ = /,. = V] |n + r) ( 



n + r . 



(31) 



Eq. f[3Tl) is equivalent to an infinite set of integral conditions for a sought for positive 
function PV^d^P) that are obtained by performing the angular integration over 6 



[z = \z\e 



ie\ 



X" 



TT 



Wrjx) 



x] 



dx = Pr{n), n = 0,1, 



X = \z\ 



(32) 



This is the Stieltjes moment problem [301 [31] ^oi Wr{x) = '7iWr{x)/JVr{x). The positivity 
requirement for Wr{x) can be always satisfied as the moment Prin), see Eq. ([22|) . is a 
product of factorials and Eq.( [32l) can be rewritten as a Mellin transform [32] of Wr{x) 
with n = s — 1 (complex s) as follows: 



x'-^Wr{x)dx = M [Wr{x)- s] 



V-1 



Wns+k) 



.k=0 



r(s + r] 



= Pr{s - 1). 

The formal solution of Eq.ljMl) is Wr{x) = Ai~^[pr(s — 1); a;], or equivalently 

JJr(s + A;) I r(s + r);x , x > 0. 

vfc=0 / 



WJx) = M~^ 



(33) 
(34) 

(35) 
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In the above formulae Ai and M.^^ denote the Melhn and inverse MeUin transforms 
respectively [32]. The positivity of Wr{x) follows from solving Eq.( l33|) by the Melhn 
convolution [10] which by itself conserves the positivity, see Refs. [28| [33| IM] for an 
elaboration of this property. 

We can offer, for any r, explicit representations of Wr{x) in terms of Meijer's G 
functions G{. . .) [9] which, via Eqs.f p3|) and fl35|) . are given by [9l [TO] : 

Wr{x) = G{\[ ],[ ]], [ [0,0,1, l,... , r-l,r-l ,r],[]],x). (36) 

2r entries 

We use the following transparent notation for Meijer's G function borrowed from 
computer-algebra systems [35] : 

G{\[list of la parameters Oj], \iist of l^ parameters jj]], 

[[list of Ip parameters (3j\, [list of Is parameters Sj]],x) = 



M 



1 ^ 



(37) 



n;=ir(7, + ^)n;Lir(i-5,- 

Note that the parameter list in the third bracket in Eq. (!36|) has 2r+l elements. Although 
apparently none of the functions of Eq. (l36p can be represented by the known special 
functions we have at our disposal a rather complete knowledge about their characteristics 
P ES]. We present on Fig. CD the plot of W;(x) for r = 1, 2, 3 . Note that aU Wr{x) are 
singular at x = for all r = 1, 2, . . .. The solutions specified by Eq. (l36|l are not unique, 
see Appendix A. 

We conclude that the states \z)r are the coherent states in the sense of Klauder. 
They are normalized, see Eq. flT3|) . continuous with label and satisfy the resolution of 
unity with a positive weight function Wr{x). 

6. Physical and statistical properties of coherent states \z)r 

In this paragraph we shall present in detail the quantities enumerated in Eqs. fl23|) - fl29|l 
calculated for different values of r. 

The probability Pr{k,x), defined in Eq.f E^ . is presented in Figs. |2] and [3] as a 
function of x for r = 1,2,3, for k = r and k = r + l, along with the Poisson distribution 
originating from the standard CS. 

The average numbers nr{x) of photons in the state \z)r are presented in Fig. HJ 
The dependence of nr{x) for fixed x as a function of r indicates increasing deviation 
downwords from the straight line n{x) = x, characterising the standard CS. This 
behaviour clearly differs from the results describing the photon-added states considered 
in [5], compare Fig. 2 of this reference. 

The Mandel parameter QM,r{x) is presented in Fig. [S] For all r the non-Poissonian 
character of the states l^;),. is explicit and is increasing with increasing r. 
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Figure 1. Plot of the logarithm of weight functions In Wr{x) , see Eg . ((36)) . for 

r = 1,2 and 3, as a function of a; = |zp. It is seen that for x > the functions of 
Ea. p6p are positive and are singular at x = 0. 



In Fig. iniwe present the plot of the metric factor Ur{x) for r = 1,2, 3. We observe 
the increasing deviation of Ur{x) with increasing r from the flat geometry of standard 
CS characterized by u{x) = 1 [2B] . 

In Fig. [7] we present the uncertainty (AX)"^ for several values of r. For sufficiently 
large x for every r (AX)^ becomes smaller 1/2 (i.e. we observe squeezing in x), although 
in Fig. Hit is only visible for r = 1 case. In Figs. [H] and [9] we present the squeezing in the 
complex plane for X and P quadratures respectively. As a guide for eye the squeezed 
regions are colored in blue. The Heisenberg uncertainty relation AX ■ AP > 1/2 is 
evidently satisfied for all values of r. It can be already seen in Figs. |8]and[9l for r = 1, 
by reading off the appropriate values for a given x. A similar behaviour persists for 
other values of r and corresponding schemes are not reproduced here. 



7. Discussion and Conclusion 



We have initiated in this work a method of generating collective quantum states of 
the Fock space via the action of certain non-unitary analogues of conventional boson 
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P,{r,x) 




Figure 2. Plot of the probability Pr{r, x) of finding the lowest allowed state \r) in the 



CS 1^)^ for 



1,2 and 3, as a function of x = |zp, see Eq. ([23)) . The dashed line is 



the probability of finding the state |0) in the standard CS, i.e. e 



displacement operator. The prescription involves first finding eigenfunctions Er{x) of 
simple classical differential operators. Then the non-unitary analogue of displacement 
operator is obtained by putting the boson creation operator as the argument in those 
eigenfunctions Er{za'^) with complex z. This procedure allows in a systematic way 
the construction of quantum states \z)r via action _E',.(2a^)|0), see Eq.lfTTl) which is the 
lynchpin of this method. In other words we depart from the exponential as a vehicle to 
create the states and look for other, related families of functions (here of hypergeometric 
type) for that purpose. In this paper we have limited ourselves to an operator of the form 



x^d"^^^ which quite naturally produces eigenfunctions of (a^ 



,r+l 



with the eigenvalue 



z. The fact that the quantum CS so obtained almost automatically fulfil the criteria of 
Klauder, is quite satisfying and may give rise to many possible extensions. 

Since our state \z)r is reminiscent of photon added-state of Agarwal et al, it is 
natural to compare the properties of \z)r with those of states of Ref. [5]. Since their 
Fock space expansions are different (compare our Eq. flT^ with Eq.(2.6) of Ref. [5]), 
one expects that their physical and statistical properties differ too. In fact this is the 
case and it appears that there is no way to relate these two families of coherent states. 
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r = 2 



0.6- 



0.5- 



0.4- 



''^'-^'•"^ 0,3 H 



0.2 



0.1 








Figure 3. Plot of the probability Pr{r + l^x) of finding the first allowed excited state 



\r + 1) in the CS \z)r for r = 1, 2 and 3, as a function of : 



, seeEq-dlSl). All the 



curves display one maximum moving away from a; = for increasing r. The dashed 
line is the probability of finding the state |1) in the standard CS, i.e. xe~^ . 



However, there are some common features: both states are squeezed, are sub-Poissonian 
in nature and possess the resolution of unity (unique for states of Ref. [5], see Ref. [15], 
and non- unique for states considered in the present work). It is interesting to observe 
that the eigenproperties of these famihes can be obtained exactly: the relevant operator 
in our case is quite natural as it is given by the boson form of differential operator whose 
eigenfunctions we consider. In the other case the appropriate operator is quite involved 
and has been obtained some time after the states have been proposed [E 
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Figure 4. Plot of the average number nr{x) of photons m the state |z)r, see Ea.p4 
for r = 1, 2 and 3, as a function of a; = |zp. 



Appendix A. 

A special feature of the solutions of the Stieltjes moment problem, Eq. flS^ is the fact that 
the Wr{x) defined by Eq. (!35|) and ( 136|) and visualized in Fig. [H are not unique solutions 
for all r = 1, 2, . . .. It means that there exist other positive solutions Vr{x) ^ Wr{x) such 
that J x^Vr{x)dx = Pr{n), n = 0, 1, . . .. In order to prove the non-unique character of 



solutions of Eq. fl32|) we apply the following sufficient condition for non-uniqueness of the 
Stieltjes moment problem; it involves both the moments Pr{n) and the solution Wr{x) 

m-. If 



Sr = EZolPrin)]-'^'"' < oo, and 



(b) the function ipr{x) 



In 



Wr{exp{x)) 



is convex for x > 0, 



then Wr{x) in non-unique [3l [371 |38l [39] • 

Concerning (a): we establish, via logarithmic test of convergence [9], that Sr < oo for 
all r = 1, 2, . . . . For (b) it suffices to prove that for fixed r the condition [■i/'rl^;)]" > 
is fulffiled. In the general case it is difficult to show it analytically, but it can be done 
relatively easy with the help of computer algebra systems [35] . 
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Figure 5. Plot of the Mandel parameter Qm.t^^) for r = 1,2 and 3, as a function 
of a; = |zp, see Eq. (p5| . The states \z)r are sub-Poissonian, because QM,r < for all 
X > 0. Since for fixed x, |(5j\/,r+i| > |QA/.r| the sub-Poissonian character is increasing 
with increasing r. 



The actual construction of such non- unique solutions, i.e. finding functions Vr{x) 
referred to above, is still an open and challenging problem, see [3l] and reference therein. 
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Figure 8. Plot of (AX)^ = X'^ — X"^ for |z)i as a function oi z = x + iy. 
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Figure 9. Plot of (AP)^ = P^ — P^ for |z)i as a function oi z = x + iy. 



